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It is shown numerically that the Coulomb self-energy of an isolated, color non-singlet source 
diverges in an infinite volume. This is in accord with the Gribov Horizon scenario of con- 
finement advocated by Gribov and Zwanziger. It is also shown that this divergence can be 
attributed to the presence of center vortices in thermalized lattice configurations. 



The energy of an isolated color non-singlet source, in a confining theory, must be infinite in 
an infinite volume, even when the usual ultraviolet divergence is regulated. We would like to 
understand this fact a little better. Color confinement is usually attributed to center vortices 
or abelian monopoles (it is now understood that these objects are related 1 ), but there is also 
another approach to confinement — the Gribov Horizon scenario, formulated in Coulomb gauge 
— which has been advocated by Gribov and Zwanziger. 2 In this talk I would like to briefly explain 
the Gribov Horizon scenario, show numerically that it accounts for the infrared divergent self- 
energy of color sources, and show also that this success is associated with the presence of center 
vortices in gauge field configurations. A more detailed exposition, on which this talk is based, 
can be found in a recent article. 3 

As a warm-up exercise, let us consider a familiar abelian example: the electrostatic energy 
of a static electric point charge, in a cubic volume of length L, with appropriate boundary 
conditions. In Coulomb gauge the Hamiltonian has the form 
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where £ is the Coulomb self-energy of a charge at point x. M is the Faddeev-Popov (FP) 
operator of the abelian theory. Let A e ^ = — d^O be a gauge transformation of a configuration 
A u fixed to Coulomb gauge. Then 
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The eigenstates of the FP eigenvalue equation M(f>^ = X n (j)^ are of course just the plane wave 
states, with eigenvalues equal to squared momenta. In a finite volume these states are discrete, 
with a lattice regularization their number is finite, and we can express the Green's function 
corresponding to M as 
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After some simple manipulations, one finds that 
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where F(X n ) = (<j)( n '\(— V 2 )|</>( ra )), and L 3 is the number of lattice sites in the cubic volume. Let 
p(X) denote the normalized density of eigenvalues. Then at large volumes we can approximate 
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the sum over eigenstates by an integral 



£ = ^j dX e^p» (5 ) 

In QED, its easy to show that p(A) = \^X/(4tt 2 ) , F(X) = X. The eigenvalues span a finite 
range, with \ m in ~ 1/L 2 and X max ~ l/o 2 , where L is the extension of the lattice, and a is the 
lattice spacing. Putting it all together, we find that that while £ has an ultraviolet divergence in 
the continuum a — ► limit, it is finite in the infinite volume L — ► oo limit. However, finiteness 
at infinite volume clearly depends on the small A behavior of p(X)F(X). If instead we had 
lim.\^ p(X)F(X)/X > 0, then the Coulomb energy would be divergent in the infinite volume 
limit. 




Figure 1: Data for full, unprojected configurations. 

In non-abelian theories, there are many gauge copies — Gribov copies — that satisfy the 
Coulomb gauge condition. The Gribov Region is the space of all Gribov copies with positive 
FP eigenvalues. The boundary of the Gribov region is the Gribov Horizon. Configurations lying 
on the Gribov Horizon have at least one FP eigenvalue A = 0. Typical Coulomb- gauge lattice 
configurations are expected to approach the Gribov horizon, in the infinite-volume limit, due 
to entropy considerations. 4 But what counts for confinement is the density of eigenvalues p(X) 
near A = , and the "smoothness" of these near- zero eigenvalues, as measured by F(X). In 
non-abelian theories the FP operator is gauge-field dependent, i.e. M = — V • D where is the 
covariant derivative, and the self-energy of an isolated static charge in group representation r is 
proportional to the quadratic Casimir C r , and to £ in eq. (5). So our approach is to calculate 
p(X) and F(X) by lattice Monte Carlo simulations of lattice SU(2) gauge theory, and extrapolate 
the results to infinite volume. 

The results are shown in Fig. 1, obtained at (3 = 2.1 for a variety of lattice sizes. We have 
applied a scaling analysis derived from random matrix theory, 5 based on the scaling of the low- 
lying eigenvalue distributions with lattice size L, to estimate that in the infinite volume limit 
we have p(X) ~ A ' 25 , F(X) ~ A ' 4 at small A. Substituting these power behaviors into eq. (5), 
we find that £ has a divergence in the infrared L — ► oo, X m i n — > limit, in addition to the usual 
ultraviolet divergence in the continuum a — > limit. In other words, the Coulomb self-energy of 
an isolated color charge is infrared infinite, by the mechanism envisaged in the Gribov Horizon 
scenario. 

We will now display a connection with the center vortex confinement mechanism. Recall 
that center vortices are surfacelike objects in the (D=4) SU(N) gauge theory vacuum, which 
can be topologically linked to closed loops. Creation of a center vortex linked to a Wilson 
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(a) Density p(A) of low-lying eigenvalues. 



(b) F(X) for low-lying eigenstates. 



Figure 2: Data for vortex-only configurations. 



loop multiplies the loop by a center element of the gauge group. The center vortex theory of 
confinement holds that the area law falloff of Wilson loops is due to vacuum fluctuations in the 
number of vortices linking the loop. 

In 1997, methods were devised for locating center vortices in lattice configurations, and, a 
little later, for removing them. This was followed by many investigations of the vortex theory in 
the lattice community, with results discussed in two recent reviews. 6 It was found that (i) center 
vortices, by themselves, account for the bulk of asymptotic string tension; (ii) the density of 
these vortices scales according to asymptotic freedom; and (iii) removing center vortices removes 
the string tension, removes chiral symmetry breaking, and sends the topological charge to zero. 

We have used the standard technique of maximal center gauge fixing plus center projection 
to separate each Monte Carlo generated configuration into two components: the vortex-only (or 
"center-projected") component, containing only the identified center vortices, and the vortex- 
removed component, in which those same vortices have been removed from the original lattice 
configuration. Each component is then transformed to minimal Coulomb gauge. Our data for 
the vortex-only configurations is shown in Fig. 2. This time a finite-volume scaling analysis 
indicates that p(X) ~ A 0±0 - 05 , F(X) « 1 as A -► 0, which again implies (from eq. (5)) an infrared 
divergence of the Coulomb energy, resulting from the vortex configurations alone. 

By contrast, when vortices are removed from thermalized lattice configurations, there is a 
dramatic change in the eigenvalue spectrum, illustrated in Fig. 3 for the 20 4 lattice volume. 
Inspection of this data reveals that the number of eigenvalues in each peak of p(A), and each 
band of F(X), matches the degeneracy of eigenvalues of —V 2 , the zeroth-order Faddeev-Popov 
operator, at the given lattice size. We know that p(X) for the —V 2 operator at finite volume 
is just a series of r5-function peaks. In the vortex-removed configurations, these peaks broaden 
to finite width, but the qualitative features of p(X), F(X) at zeroth order, i.e. the absence of 
confinement, remains. 

These numerical results establish that the Coulomb self-energy of a color non-singlet state is 
infrared divergent, due to the enhanced density p(X) of Faddeev-Popov eigenvalues near A = 0. 
This supports the Gribov-Zwanziger picture of confinement. It also appears that the confining 
property of the FP eigenvalue density can be entirely attributed to center vortices, since (i) 
enhancement of p(X), F(X) is found in the vortex-only content of lattice gauge configurations; 
while (ii) the confining properties of p(X), F(X) disappear when vortices are removed from the 
lattice. 

We conclude with two further facts about center vortices and the Gribov horizon, stated here 
without proof: First, vortex-only configurations have non-trivial Faddeev-Popov zero modes, 7 
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(a) Density p(A) of low-lying eigenvalues. 



(b) F(X) for low-lying eigenstates. 



Figure 3: Data for vortex-removed configurations. 



and therefore lie precisely on the Gribov horizon, which is a convex manifold in the space of 
gauge fields, both in the continuum and on the lattice. Secondly, vortex-only configurations 
are conical singularities on the Gribov horizon. 3 It appears that center vortices have a special 
geometrical status in Coulomb gauge, although the physical implications of this fact are not yet 
fully understood. 
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